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Abstract 

For one dimensional SU(n) Hubbard model, a pair of Lax operators are derived, which 
give a set of fundamental equations for the quantum inverse scattering method under 
both periodic and open boundary conditions. This provides another proof of the inte- 
grability of the model under periodic boundary condition. 

PACS Numbers: 46.10. +z, 05.40.+j, 05.60. +w 



e-mail: yue@phys.ocha.ac.jp 



Integrable strongly correlated electron systems have been an important research subject 
in condensed matter physics and mathematics. One of the significant models is the 1-D 
Hubbard. The exact solution was given by Lieb and Wu Jl|. However, the integrability was 
shown twenty years later by Shastry, Olmedilla and Wadati 0, [3|] . The integrability and the 
exact solution of the system under the open boundary condition were discussed by several 
authors f|, ||. The Lax pair was first given by Wadati, Olmedilla and Akutsu ||. Recently, 
Maassarani and Mathieu have constructed the hamiltonian SU(n) XX model and proved 
its integrability [[7|. Considering two coupled SU(n) XX models, Maassarani succeeded in 
generalizing Shastry's method to SU(n) Hubbard model ||. Further, he solved the Yang- 
Baxter equation to prove the integrability of one dimensional SU(n) Hubbard model |J. 
( It is also proved by Martins for n = 3,4 |1(J.) In this paper, we apply the quantum 
inverse scattering method to 1-dimensional SU(n) Hubbard model and derive the explicit 
form of the Lax pair, which gives another proof of the integrability. It is worthy to note 
that straightforward application of the reflection matrix method to study the problem of the 
integrable open boundary in the present SU(n) Hubbard model would encounter difficulties 
due to the fact that the R matrix in reference [EJ does not satisfy the crossing symmetry 
condition and the lack of invertibility of a certain matrix. 

The Lax pair formalism will give an effective method for such a system. The 1-dimensional 
SU(n) model in the Schrodinger picture is given by 

N n-1 t j 2 N 

U \ ^ \ ^ ( rpna rpan , rpan rpna , rpna rpan , rpan rpna \ , \ " r~i r~i /-i \ 

ri= 2^ [■ b <r,k- b <r,k+l + Ar.fc-^o-.fc+l + ^r.fc-^r.fc+l + ^r.fe-^r.fe+l ) + 2^ C ^T,k, (1) 

k=la=l 4 k=l 

where E^ k {a = a, r) is a matrix with zeros everywhere except for a 1 at the intersection of 
row a and column f3: 

(E a %m = «• (2) 

The subscripts a, k stand for two different E operators at site k (k = 1,...,N). The 
n x n diagonal matrix C is defined by C = J2 a <nE aa — E nn . The Hamiltonian enjoys the 
(su(n — 1) © w(l)) CT © (su(n — 1) © u(l)) T symmetry. The generators are 

N N 

J f = 2 E a% and K a = Yl C ^k, a,(3=l,...,n-l, a = a, T. 

k=l k=l 

In the rest of this letter we discuss various operators in the Heisenberg picture. When we 
deal with the operators corresponding to the matrices E and C in the Schrodinger picture, 
they are denoted by adding a hat ( " ) to the corresponding matrix: 

Q(t) = e mt Qe- mt , Q = E or C. 
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In the following, we do not indicate the time dependence of the operators. Applying this 
method to the Hamiltonian ([[]), we find 



dt 



% 

f3<n 



/ j ^• C/ o,fc- C/ a 1 fc+l -^a.fc-^a.fc+l ^ -^a.fc-^a.fc-l - C/ a,fe- C/ a,fe-l / 
3<n 

dP na Tin 2 

a M _ ■ f fpfia fpnfi , fpfia fhnf3 \ ■ ( fpnn fpna , fpnn fpna \ ■ ulb fpna/y 

i. — 2 V^a.fc-^a.fc+l + ^a,k^a,k-l J ~ 1 {^a,k^a,k+l + ^a.fc-^a.fc-l J - ! ~ 

/3<n 

( -//? cm Tin 2 

<hk _ ■ ( fpnn fpan , Trmn pan \ • / p>a/3 77i/3n _j_ 77ia/3 7ri/3n \ . - U n fpan A 

j t — 1 [f'a.k&afi+l + ^a,k^a,k-l J ~ 1 2^ V^a.fc-^a.fc+l "+" ^a,k-^a,k~l J "+" '"^o,^",; 



/3<n 



■ / pan pan 77in/3 , 77in/3 pan pan pn/3 \ /qN 

^ 1 y^a^a^+l ~ -^a.fc-^a.fc+l "+* • c 'a,fc- c 'a ) fc-l ~~ -^a.fc-^a.fc-l J 5 l J J 

where (a = a, r) and (a = r, a) . For an infinite system, it is not necessary to specify 
the boundary condition. However, one should understand E aQ is equal to E a<N under the 
periodic boundary condition. In an open boundary system, E a p and -E a ,iv+i in the r.h.s. 
must be regarded as vanishing. 

Let us first consider the degenerate case U = (the SU{n) XX model JFj). In this case, 
the equations of motion ([|) decouple into two identical sets of equations for a and r. Let us 
introduce the L-operator for each of them: 

L , fe (A) = cos(A) S$> + sin(A) T$ + U$> , (4) 

in which the blocks , and are defined by 

c(+) I f? a P JpP a \ i fpnn rpnn 

J a,k — V 2-~i ^a.k^a.au) ' JZ 'o,fe- C/ a,oMJ 
a,/3<n 



_ \ I fpnn rpaa _i_ paa rpnn 
a,k / j [ a,k a,au ' a,k a, at 

a<n 

rr(+) / fpna rpan , fpan rpna \ /r\ 

^a,fc - { E a,k E a,au + E a,k E a,au ) • l 5 J 

a<n 

They satisfy various identities inherited from the definition in terms of E af3 's (||): 

_ rp{+) q (+) _ q(+) t A+) _ rr(+) Q (+) _ n r (+)rr(+) _ rr(+)m(+) _ rr(+) / R x 
°a,k 1 a,k — 1 a,k°a,k —°a,k u a,k — u a,k°a,k ~ U 5 1 a,k u a,k ~ u a,k 1 a,k ~ u a,k ? • • • 5 • W 

We also introduce the M-operator for each species: 

M a)fc (A) = ^ {^l^fc-^fc-l + A ^ E Z^ E a%-\)\ E a n au 

0<n 

+ J2 B {(41 + + $5 + 

/3<n 



+ E DiKZEZU + DzE-E™^ + D 3 (E^E n a ^ + E^E^) E™ u 

a<n [ 0yta<n J 

+ E E ^KRTi-^fitjC (7) 
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where A, B, D, F are as yet undetermined functions of a (spectral) parameter A. The matrices 
E a ,au are the constant matrices with the same definition as E a ^, k = 1, • • • , N. The subscript 
au stands for the auxiliary space instead of the quantum space. Thus L a ^ and M a ^ are nxn 
matrices in the auxiliary space. We want to rewrite equations (||) in a matrix Lax pair form 

= M a , fc+1 (A)L a , fc (A) - L ajk (X)M a>k (X). (8) 

Substituting equations (||) and (|7J) into the above Lax-pair form (|8|), we find the solution 
A\ — Di — i + ztan(A), A 2 = D 2 = i — itan(A), 

(9) 

D 3 = i, B = — i/cos(A), F = ztan(A). 

From now on we do not denote the A-dependence of L a ^ and M a ^ for brevity. The transfer 
matrix for SU(n) XX model with N sites can be defined by 

T a ,xx = L a N ■ ■ ■ L 0j i, 

which satisfies 

dT 

,N+iT a xx — T a xxM a l . 

at 

From this it is standard to show that the trace of T a ^xx is independent of time under periodic 
boundary condition. So the SU(n) XX model is integrable. This Lax-pair is very important 
for the construction of the Lax-pair of SU(n) Hubbard model. Notice that its integrability 
was first given by Maassarani and Mathieu in the framework of Yang-Baxter relation 171]. 

Now, let us consider U ^ case. In terms of the above L a ^ and M a> fc, we can rewrite 
equations © as 

dL a j Un 2 

— M a)j+1 L a j - L a jM a j + i—^-[L a j, C ajau Ca,j), (10) 

in which the last term in the r.h.s. manifests the coupling between the two species a and r. 
Using the relation [L a ,j, C aiCm ] = -[L a j, C a>j ], we obtain 

dL a ■ Un 2 

^ — M a j +1 L at j — L a jM a j — i—^-[L a j, C a jCa,j}- (11) 

We define the following operators for the coupled system 

Lj = L ad L T j , Mj = M ad + M T j , (12) 

then equation ( [TTl) can be written as 

dT ~ ~ ~ ~ Un 2 ~ 

= M J+1 L 3 - L 3 M 3 - i—[L j} (■„,(■;,.. (13) 



Now, we want to rewrite the last term in the above equation so that we could obtain the 
Lax-pair form for the coupled system in a similar form to (|8]). Following the method given by 
Wadati et al for SU{2) Hubbard model, we introduce a "rotation" matrix by the w(l) 
charge: 

h 

I au = cosh(/i/2) + smh(h/2)C^ au C T , au = exp{-C CTiau CV, a «}, (14) 
where h is a free parameter to be determined later and the "rotated" operators are 

Cj = IauLjlau, Mj = IauMjlau- (15) 

By this the fundamental blocks , and Uj^ are mapped to their 'anti-symmetric' 
counterparts S { ~ k } , T o y and U$: 

q(~) C< c(+) q(+)/^ / \ ^ f? a @ J?P a \ fpnn rpnn 

°a,k — ^ 'a,awJ a ,k ~ J a,k u 'a,au ~ [. ^a^a&u) ^tt.i^aui 

rp\—) ft rpvr) rp{+) f~t / fpnn rpaa haa rpnn \ 

1 a,k ~ ^a,auJ- a ,k ~ 1 a,k ^"a,au — 2_^i \ I1 'a,k I1 'a,au ■ Cj a,k I1 'a,auJ > 

a<n 

rr(-) _ n r/(+) _ _Tj( + )n - ( fpnarpan _ fpan rpna \ (if>\ 

u a,k ~ ^a.au^a.k ~ u a ,k Ly a,au — 2—i \ a,k ^a^au 111 a,k 11 ' a,au ) ' \ lu ) 



a<n 



which also satisfy identities similar to those given in (^). Then equation ([03) becomes 
dC ■ Un 2 

= ilMj+jJc, - c 3 m 3 - i— [£,, r ;./' ,;• (17) 

Using the definition of I au and /A, we obtain 

llM j+1 lJ = M j+1 + Q j+1 + Qj, (18) 
Qj = --^^smh(h)(ul-)C T , au + U^C a>au ), (19) 



cos(A) 



and 



dt 



(Mj+x + Qj+i^j-CAMj + Qj) 
Un 2 

+CjQj + QjCj — i [Cj, C a jC T j]. 
Detailed calculation shows the last line in the above equation to be 



(20) 



2 



Un 



I 



1 \ . , , , s . Un 2 a + b' 



a{ b -—b) Sillh(2A) +i —— b \ [C " C °"' C -"" ] (21) 



.Un 2 2ab 



L auj 
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where a = cos(A), b = sin(A). In deriving equation (^), use has been made of the following 
identities 



(a — b)[Cj, C a> jC Ti j] + (a + b)[Cj, C a ^ au C T ^ av \ 
= -4abl au (U ( J (Stf + T#) + U# (Slj> + Ttfj) I au , (22) 

huQjlau = -^stoh(2A)(E7#^ 

which are straightforward consequences of the identities (§) . In expression (|2"ID let us choose 
the parameter h by 

Un 2 Un 2 

sinh(2/i) = —2ab = — sin(2A), (23) 
so that the second line in the r.h.s. of ( pij ) vanish. Then the Lax-pair can be written as 

^ = isy,r„ (24) 

Un 2 i f 1 \ a + 6l 
B j = M j + Q j +i—l2b(b-—)+^\C a>au C Tjau . (25) 

Here, A4j,Qj are defined by equations fllSD and (^9|) , respectively. It is clear that the 
condition (p3|) guarantees the parametrizability for SU (n) Hubbard model, which was first 
introduced in || through the Yang-Baxter equation. As we know, the standard method to 



construct the open system is to study the reflection equations ||1 1|| . The R matrix must enjoy 
unitarity and crossing symmetry. For SU(ri) Hubbard model, the R matrix in does not 
have the crossing symmetry and the invertibility of a related matrix is lacking. In the Lax 
pair formalism, however, the detailed properties of the R matrix do not come in. This is 
another way to study the integrability of the open boundary systems. Therefore, the Lax- 
pair derived in the present letter would be a useful tool for analyzing open systems. We will 
consider its applications elsewhere. 
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